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ABSTRACT 

We investigate the orbital dynamics of a barred-spiral model when the system is rotating 
slowly and corotation is located beyond the end of the spiral arms. In the characteristic of the 
central family of periodic orbits we find a “bistable region”. In the response model we observe 
a ring surrounding the bar and spiral arms starting tangential to the ring. This is a morphology 
resembling barred-spiral systems with inner rings. However, the dynamics associated with this 
structure in the case we study is different from that of a typical bar ending close to corotation. 
The ring of our model is round, or rather elongated perpendicular to the bar. It is associated 
with a folding (an “S” shaped feature) of the characteristic of the central family, which is 
typical in bistable bifurcations. Along the “S” part of the characteristic we have a change 
in the orientation of the periodic orbits from a xl-type to a x2-type morphology. The orbits 
populated in the response model change rather abruptly their orientation when reaching the 
lowest energy of the “S”. The spirals of the model follow a standard “precessing ellipses flow” 
and the orbits building them have energies beyond the “S” region. The bar is structured mainly 
by sticky orbits from regions around the stability islands of the central family. This leads to 
the appearance of X-features in the bars on the galactic plane. Such a bar morphology appears 
in the unsharp-masked images of some moderately inclined galaxies. 
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1 INTRODUCTION 

In dynamical systems a “bistability situation” usually refers to 
cases where a system has two stable equilibrium states. In a bi¬ 
furcation diagram the curve of steady state displays an “S” shape 
as a certain parameter of the system varies. The “S” is delimited 
by two saddle-node bifurcation points. Between them we have two 
stable and one unstable ste ady states (see e.g. lAngeli et"a3l2003l : 
lLvnchl2007l : IStrogatzl2014l) . The corresponding situation in Hamil¬ 
tonian Galactic Dynamics is depicted in the characteristic of a fam¬ 
ily of periodic orbits as two successive tangent bifurcations (see 
e.g. IContopoulosf[2004 ) facing opposite directions. These two bi¬ 
furcations share the same unstable branch. In other words the char¬ 
acteristic folds twice as the varying parameter, usually the Jacobi 
constant E j , i ncreases. Fold i ngs of the c haracteristic have been en¬ 
countered by ISkokos etafl d2002ah and ISkokos etafl d2002bh . in 
3D Ferrers bar potentials. However, the way they affect the face- 
on morphology of a model has not been examined in those papers. 
Nevertheless, it was clear that the foldings of the characteristics 
affect to a larger degree slowly rotating models. In this paper we 
present the implications of the presence of such a folding of the 
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characteristic of the main family of periodic orbits for the dynam¬ 
ics of a slowly rotating barred-spiral potential. 

Slowly rotating models of disc galaxies have been proposed 
in the past to describe the dynamics of normal (non-barred) spi¬ 
ral galaxies with open spiral arms. In these models (stellar and 
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Contrarily, in barred galaxy models_ corotation is usually 
placed close to t he end of the bar dContopoulosi ri980). Recently 
iFont et all d2014h presented a list with 32 barred galaxies in which 
they estimated the ratio of the corotation to the bar radius, RJRb , 
to be between 0.94 + 0.08 < RJRb < 2.1 + 0.5. Model bars 
ending well inside corotation have been found in N-body simu¬ 
lations JCombes & Elmegreenlll993l) as well as in response mod- 
els of barred potential s derived from near-infrared observations 
dRautiainen et al.lE008l) . In all these studies, slowly rotating bars 
are associated with late-type barred-spiral galaxies. It is generally 
believed that bars in barred galaxies are supported by the xl family 
of elliptical, stable, p eriodic orbits, which extends^along the ma¬ 
jor axis of the bar iContonoulos & Papavannopouloslll980h . or, in 
the three-dimensional case, by the corresponding families of the 
xl-tree, i.e. by xl together with 3D families bifurcated at the ver- 
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tical resonances, (ISkokos et al.112002a!) . The bar is built by trapping 
quasi-periodic orbits along the orbits of these families. Deviations 
from this orbital behaviour have been proposed, pointing to bars in 
which other families than xl play a significant and perhaps a lead¬ 
ing role in the building of the bars. Such behaviours are favoured ei¬ 
ther in slo wly rotating bar models in which the bars end well befor e 
corotation dPetrou &Papavannopoulosfl986l : ISkokos et ali 2002b). 
or in b ars with large major to minor axis ratios (Kaufma nn & Patsisl 
l2005h . |Pasha & Polvachenkol i 19941) even claim, that in slowly ro¬ 
tating bars of the type described bv lLvnden-Belll dl979h . there is a 
better matching of the outer-to-inner ring radial ratio, than in stan¬ 
dard fast rotating bars. Finally. IPatsis et al.l ( 201 0) have shown that 
bars of ansae-type can be supported mainly by chaotic orbits in a 
(fast rotating) model based on the potential of NGC 1300, estimated 
from near-infrared observations. 

The implications of slow rotation for the orbital dynamics in a 
barred-spiral model, i.e. when the spiral component is explicitly in¬ 
cluded in the potential, h ave not been extensively studied. Recently 
ITsigaridi & Patsisl d2013h have presented a barred-spiral model, ro¬ 
tating with a single pattern speed, characterized by a ratio RJRb , 
about 2.9. This case, with Q p = 15 km s -1 kpc -1 , has been consid¬ 
ered as “general” since several dynamical mechanisms cooperated 
in forming the obtained barred-spiral response morphology. The ac¬ 
tion of two different dynamical mechanisms led to the formation of 
an inner barred-spiral structure surrounded by an oval-shaped disc 
and an outer set of spiral arms beyond corotation. However, if we 
decrease further the pattern speed, there are even more considerable 
changes in the orbital dynamics of the system. The pattern speed, 
Q p , is the most important parameter for the resulting response mor¬ 
phology of the model. 

Besides the case presented in Tsigaridi & Patsisl d2013h . we 
have studied a series of models with even smaller Q p down to 
10 km s -1 kpc -1 . As the pattern speed decreases from Q p = 
15 km s -1 kpc -1 , we have a rather new orbital behaviour, which 
shapes a different barred-spiral response with distinct morphologi¬ 
cal features. For example, despite the fact that the bars in all studied 
models have comparable sizes, as Q p decreases the orbital dynam¬ 
ics of the central fami ly of per iodic orbits changes from that of the 
typical case iContopoulos & Grosbcfll 119891) . In the present paper 
we present these changes, we describe the morphological features 
that are encountered in extremely slowly rotating models, and we 
discuss its relevance to some morphological features appearing in 
some images of barred-spiral galaxies._ 

Our general model is the same as in lTsigaridi & Patsisl d2013h . 
It is based on a modified version of a potential estimated from 
the near-inf rared photometry of the late type barred-spiral galaxy 


NGC 3359 teoonvasaidl2003 : IPatsis et alljoo^l . We used it 
general barred-spiral potential with Q D being a free parameter. It is 
worth noticin g tha t Elmegre cn & Elmegreen Il985 ), based on ob¬ 
servation bv lGottesmanlil982l) . estimate the ratio of the semi-major 
axis of the bar to the length of the rising part of the rotation curve 
for NGC 3359 to be Rbar/Rrise = 0-75. According to this work, the 
fact that the bar ends before the velocity curve turns over, is an in¬ 
dication for slow bar rotation. Despite the difficulties in estimating 
the inner parts of the rotation curves of barred galaxies, this prop¬ 
erty is in general associated with slowly rotating bars, like the bars 
of the models we study. In the s eries of our response mo dels, since 
Rrise « 4.6 kpc (see figure 3 in Tsigaridi & Patsis 2013) this ratio 
varies between 0.619 and 0.641. In the specific model of the present 
paper this ratio is about 0.62, close to the lowest value. This means 
that our potential is adequate for studying slowly rotating models 
in general. 
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Figure 1. The stellar response model with Q p = 11.5 km s -1 kpc -1 . The 
system has completed 10 pattern rotations, rotating counter-clockwise. The 
bar extends roughly along the y-axis. Corotation is at R c ~ 10 kpc 


The structure of the paper is the following: In Section |2] 
we briefly present again our potential (for more details see 
[Tsigaridi & Patsisll2013l) . The results of our study are described in 
Section [3] and refer to the building of the response features, which 
are the ring, the bar and the spirals. These results are discussed in 
Section [4] and in Section [5] we enumerate our conclusions. 


2 SUMMARY OF MODEL PROPERTIES 

The m odel has been extensively described in ITsigaridi & Patsisl 
J2013I) . It is a two dimensional model of the general form 

$</, (p) = d>o(r) + ®mc(r) cos (mtp) + 9 ms (r) sin ( nup ) (1) 

m= 2,4,6 

The components O 0 , ® mc , and of the equation above are given 
as polynomials of the form n = 0... 8. The radial variation 

of the perturbation force normalized over the radial axisymmetric 
one is given in figure 1 in lTsigaridi & Patsisl d2013h . 

The equations of motion are derived from the Hamiltonian 

H= l -{x 2 + y 2 ) + ®(x,y) - ^(x 2 + v 2 ) = Ej, (2) 

where ( x , y) are the coordinates in a Cartesian frame of reference 
rotating with angular velocity £l p . ®(x,y) is the potential m in 
Cartesian coordinates with the bar aligned approximately with the 
y-axis, Ej is the numerical value of the Jacobi constant, hereafter 
called the energy, and dots denote time derivatives. 


3 SLOWLY ROTATING MODELS 

By varying Q p between 10 < Q p < 30 km s _1 kpc -1 , we obtained 
always a kind of barred-spiral response. In this range of pattern 
speeds the corotation radius of the models, R c , takes values between 
12 > R c > 4.3 kpc respectively. Nevertheless, while the pitch angle 
of the response spirals varied considerably in models with different 
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Figure 2. The characteristics of xl' and its bifurcations for Q p = 
11.5 km s -1 kpc -1 . The dashed line, labelled “ZVC”, is the curve of zero ve¬ 
locity. Continuous black parts of the characteristic indicate stability, while 
dotted parts, red in the online version, instability. The E b] line is indicated 
with an arrow in the upper right corner of the figure. 


pattern speeds (it was larger for lower pattern speeds), the radius 
of the response bar varied only between 2.85 < R b < 2.95 kpc. 
For Q p > 30 km s -1 kpc -1 the size of the response bar was clearly 
decreasing. For example, for Q p = 35 km s -1 kpc -1 we estimated it 
to be about 2.45 kpc. 

The changes that are introduced in the dynamics of the sys¬ 
tem as Q p decreases are reflected in ch anges observed in the (E j,x) 
characteristic curve (for a definition see llContopouloslEool . section 
2.4.3) of the central family. In typical cases of barred galaxy mod¬ 
els, the xo initial condition of the central, xl, family is increasing 
with Ey between the inner Lindblad Resonance (ILR) and the 4:1 
resonance dContopoulos & Grosb0llll98H Il989h . This is the case 
also for the present model for Q p > 16 km s -1 kpc -1 . However, 
for slower rotating models, i.e. for Q p <; 15 km s -1 kpc -1 , we ob¬ 
serve a folding of the characteristic curve well before corotation in 
the (E j,x) diagram. The AEy range over which we have the folding 
in a model increases with deceasing Q p . We call this feature, the 
“S”. As we will see, foldings of the characteristics introduce in the 
system new orbital dynamics accompanied by new morphological 
features. For this reason we call just in this work the central family 
of our model xl' in order to distinguish it from the standard xl fam¬ 
ily, which has as members only elliptical periodic orbits elongated 
along the bar. 

In the present paper we describe these changes in a typical 
case with Q p = 11.5 km s -1 kpc -1 (chosen so that we have R c « 
10 kpc). For this Q p the morphological features associated with 
slow rotation dominate and this facilitates the description of the 
relevant dynamical mechanisms. 

The stellar response model for the potential 0 with Q p = 
11.5 km s -1 kpc -1 is given in Fig. [l] The response bar length is 
R b « 2.85 kpc. The system has completed 10 pattern rotations. 
Initially the particles have been distributed randomly within a disk 
of radius r max -\\ kpc with velocities securing circular motion in 
the axisymmetric part of the potential, O 0 (r). At the beginning of 
the simulation the amplitudes ® mc and ® mv of the perturbing term 
grow linearly with time from 0 to their maximum value within two 
system periods (“time dependent phase”) and after that they remain 
constant. The time-dependency of the amplitudes in the beginning 
of the simulation secures a smooth response of the system, since 
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Figure 3. Periodic orbits along the characteristic of xl' for Ey = -55000 (a), 
-45000 (b), -37000 (c), -33000 (d),(e),(f) (the three xl' representatives at 
the “S” region), -30000 (g) and -27000 (h). Note the different scales of the 
axes. 


the initial velocities are for circular motion in the axisymmetric 
term ® 0 (r). 

The snapshot has been converted to an image using the ESO- 
Midas software. In the central part we observe a rectangular shaped 
bar. The contrast of the image has been chosen such as to allow us 
to clearly see that inside the rectangular shaped bar appears an “X” 
feature. The bar, with the embedded in it X-feature, ends on an oval 
(pseudo)ring structure. This ring has a certain width. It is almost 
round at its inner boundary, which is attached to the bar, while the 
ellipticity of the orbits that form it increases outwards and becomes 
maximum at its outer part, which coincides with the beginning of 
the spiral structure. At this point the ring is elongated in a direction 
roughly vertical to the bar (at an angle about 15° with respect to 
the x-axis). The length of the semimajor axis of this ring is about 
3.6 kpc. On the sides of the bar, inside the ring, the response surface 
density is very low. The double armed open spiral structure starts 
tangentially from the ring having a pitch angle i p « 35°in the outer 
parts (it is not logarithmic). The spiral arms have a sharp bifurcation 
at r « 7.8 kpc. 
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Figure 4. The evolution of the surfaces of section in the ”S” region. In (a) we focus in the ”S” region of the characteristic. Continuous black parts of the curve 
indicate stability, while dotted (red in the online version) indicate instability. The vertical line at Ey = -34000 denotes the energy at which we give the (x, x) 
surface of section in (b). With increasing energy between A, A' and B, B' in (a) the size of the left stability island decreases, while the size of the right one 
increases. The arrows point to the orbits depicted in Fig. [5] 


The (E/,x) characteristic curve of xl' and its bifurcations, as 
we approach L 4 , is given in Fig. 0] Continuous black parts of the 
curves indicate stability, while dotted (red in the online version) 
indicate instability. The stable Lagrangian points L 4 and L 5 are at 
energy E L45 = -22196 (Ey is always given in units of km 2 s -2 ). 
L 4 in Fig. EJis at the local maximum of the curve of zero velocity 
at the right side of the figure. For the unstable Lagrangian points 
Li and L 2 we have E Ll 2 = -23360. At this energy we have drawn 
a solid vertical line in Fig. 0| and we indicate it with an arrow la¬ 
belled “E Ll ”. The folding of the characteristic curve, which we call 
“S”, occurs between -34500 < Ey < -31800 (Ey ~ -34500 is the 
energy at point “A” and Ey « -31800 is the energy at point “B”). 

Since there is no gap or discontinuity in the curve, we con¬ 
sider the whole curve belonging to a single family of periodic or¬ 
bits, namely the xl' . The xl' family is stable close to the centre of 
the model, then it becomes unstable for -57000 < Ey < -50436 
(with a small stability interval -55000 < Ey < -54336). After 
that it remains practically stable until the region of the “S”. For Ey 
> -31800 (Ey at point “B”), which we consider as the end of “S”, 
the upper branch of the characteristic remains stable until the en¬ 
ergy Ey « -27000. Beyond this energy follows a tree of bifurcating 
families. 

Up to this point, the morphological evolution of the periodic 
orbits along the xl' characteristic curve shows an interesting vari¬ 
ation. In Fig. [3] we give successively characteristic periodic orbits 
as the energy increases. For Ey ^ -54000 the xl' orbits are typ¬ 
ical xl ellipses, elongated along the y-axis, as in Fig. [3k- Then 
for -54000 ^ Ey ^ -38000 they develop loops at their apocentra 
(Fig. 0k). As energy increases further, for Ey > -38000, the loops 
at the apocentra vanish again (Fig. 0k) and the morphology of the 
periodic orbits as we approach the “S” region is as in Fig. 0]i In 
the “S” region we have three simple periodic orbits at each energy. 
Two stable and one unstable. At the local maximum of “S” (Ey 
« -31800), where the characteristic turns to the left (point “B” in 


Fig.0, the orientation of the elliptical periodic orbits starts chang¬ 
ing, while, simultaneously, they become unstable. Their major axes 
lean more and more towards the x-axis as energy decreases. For 
example the orbit in Fig. [3f, at Ey « -34500 (point “A ”), the ma¬ 
jor axis of the periodic orbits is close to the x-axis (minor axis) 
of the system. In other words, moving along the unstable segment 
of the characteristic from “B” to “A” w e chan ge from a xl- to a 
x2-like (see e.g. lContopoulos & Grosb0lill[l989 ) orientation. Then, 
moving again from Ey ^ -34500 to the right along the upper stable 
branch of “S”, the orbits become stable and of “x2-type” until Ey 
« -30000 (Fig.0k). Beyond this point we have periodic orbits like 
the one in Fig. 0k- If we plot together successive xl' periodic or¬ 
bits with increasing Ey for Ey > -28000 we observe that their major 
axes start tilting towards the y-axis this time, building a “precessing 
ellipses pattern” |Patsisll2009l ) that can be considered as the back¬ 
bone of a spiral structure extending to larger distances. The spiral 
structure is discussed below in section [3731 


3.1 The ring 

In Fig. 0| we focus into the “S” region. The corresponding part of 
the characteristic is given in Fig. 0k and is included between the 
energies of‘A”, ‘A'” and “B”, “B'\ “A”, “A'” are at Ey « -34500 
and “B”, “B'” are at Ey « -31800. At any energy between, e.g. for 
Ey — -34000, where we have drawn a vertical line in Fig. 0k, we 
have three periodic orbits. There is always an unstable xl' periodic 
orbit with an intermediate inclination between the stable xl-like 
(lower part of the characteristic) and the stable x2-like (upper part 
of the characteristic). The situation can be described as a case of a 
bistable bifurcation. We can say that at points “A” and “B” we have 
a direct and an inverse tangent bifurcation. 

In Fig. 0k we give the (x, x) surface of section for the energy 
Ey — -34000. It is a typical case of a surface of section in the “S” 
region. At the centre of the left stability island we have the “xl” 
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Figure 5. (a) Two characteristic quasi-periodic orbits from the ”S” region of xl' supporting a xl and a x2 flow. Both are at Ej =-34000. They are overplotted 
on our response model. The location of their initial conditions on the (x, x) surface of section is indicated in Fig. (b) Three orbits from the chaotic region 
between the stability islands of Fig. m ? plotted in the central region of the model. Instead of a ring they contribute to the formation of a bulge-like central 
component. 



Ej 



Figure 6. (a) A histogram showing the E j distribution of the particles with 
2.5 < r < 3.8 kpc in Fig. [I] (b) The (x, x) surface of section at Ej =-32000, 
which is a local maximum in (a). 


periodic orbit, while at the centre of the right stability island we 
have the x2-like ellipse. The extent of this cross section, as well as 
of all other similar figures we discuss in our paper, is limited by the 
zero velocity curve, ZVC (drawn in Fig.[4k). Moving from “B” to 
“A” along the unstable part of the characteristic in Fig. EK the area 
of the right stability island in the surface of section (Fig. [4J?) de¬ 
creases, while the area of the left one increases (Fig. ij) is close to 
the left border of the “S”). At E j ^ -33000 the width of the chaotic 
zone between the two ordered regions becomes minimum. The role 
of the x2-like orbits is even more emphasised as Ej increases be¬ 
yond “B'” along the characteristic in Fig. [4j). From “B'” until E j 
« -30000, the only simple periodic orbits existing are almost per¬ 
pendicular to the bar. Since they are stable they may attract around 
them quasi-periodic orbits. If these orbits are populated, they will 
support a x2 flow. 

By inspection of Fig. [3] we can understand the association be¬ 
tween the xl' periodic orbits and the non-periodic orbits that have 
been populated in the energy interval of the “S” region (Fig. [4j 
in order to give the response morphology depicted in Fig. Q] We 
remind that the initial conditions of the particles in the response 
model are distributed randomly on a disc of radius 11 kpc and take 
velocities for circular motion in O 0 in Eq. (1). The orbits that popu¬ 
late the response model have as initial conditions the positions and 
the velocities of the test particles at the end of the “time dependent 
phase” of the model that lasts two rotational periods. In general the 
xl' periodic orbits with small initial x 0 values (0.6 ^ x 0 ^ 1.1) 
are orientated along the bar (Fig.[3^-d). They are ellipses with their 
major axes roughly aligned with the major axis of the bar. How¬ 
ever, the projections of such orbits from the “S” region, as well as 
the projections of the quasi-periodic orbits trapped around them, on 
the major axis of the bar (y-axis) are larger with respect to the bar 
of the response model, which is included inside the ring (Fig. [lj. 
If we consider quasi-periodic orbits from the left stability island of 
the surfaces of section in the “S” region they always exceed the size 
of the ring of the response model. Contrarily, the size of the quasi- 
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Figure 7. In (a) we give the Ey distribution of particles located at r < 
2.5 kpc in our response model (Fig. \T\ The mode is at Ey = -40750. In 
(b) we give the (x, x ) surface of section at the energy of the mode. 


periodic orbits from the right stability island is such as to contribute 
to the formation of the ring, without exceeding the dimensions of 
the response feature for all Ey values along the “S”. Two quasi- 
periodic orbits for Ey = -34000 that illustrate the above statement 
are given in Fig. [5b- The locations of the orbits drawn in Fig. [5b 
in the corresponding stability islands are indicated with arrows in 
Fig-lib- The one with the small x 0 in Fig. [4b is indicated with “xl” 
in Fig. [5^. 

Chaotic orbits associated with the unstable branch of the “S” 
are also not populated in the response model (Fig. |T|). In the case 
depicted in Fig. [4b, such orbits are found in the chaotic zone be¬ 
tween the two stability islands. Integrated for 8 rotational periods, 
they fill the whole central region of the model, including the rather 
empty areas between the bar and the beginning of the spiral arms, 
inside the ring. The presence of these orbits in the response model 
would build a kind of response bulge, instead of a ring. This can 
be seen in Fig. 05, where we plot on top of the model, three orbits 
from the chaotic region with initial conditions (x, x) = (2.1,-5), 
(1.7, -35) and (1.9, -35). We observe that they fill the interior of 
the ring, unlike with what happens in the model. So we can exclude 
them from the orbits contributing to the response morphology ob¬ 
served in Fig.[l] 

The identification of the ring with the x2-like orbits in the “S” 
region (upper stable branch of “S”) can be verified by considering 
the distribution of the energies of the particles located in the annu¬ 
lus with 2.5 < r < 3.8 kpc in the response model. This is roughly 


the region of the ring. The radius r = 2.5 is reaching the inner part 
of the drawn “x2” orbit in Fig. [5b, while r - 3.8 reaches the drawn 
apocentra of the “xl” orbit and the beginning of the spiral arms in 
the same figure. This energy distribution is given in Fig. [6b. We 
observe that most of the particles in the ring have energies in the 
range -35000 ^ Ey ^ -28000. Particles with Ey > -28000 are in 
quasi-periodic orbits trapped around the stable periodic orbits par¬ 
ticipating in the beginning of the “precessing ellipses” pattern that 
builds the spirals (see also Section [33] below). Their small contri¬ 
bution in the distribution of the energies in the annulus can be seen 
at the right part of Fig. [6b for Ey > -28000. For Ey > -32000 we 
do not have xl' orbits with orientation along the major axis of the 
bar. This means that if we find in the ring area particles in bar¬ 
supporting xl-like orbits they should have Ey < -32000. However, 
as we have seen, we have excluded the presence of xl-like orbits 
from the lower branch of “S” in the response model, since their size 
exceeds the size of the ring. So, according to Fig. [6b, if such orbits 
exist, they will be a minority with Ey < -34500. We also note that 
as we move in the lower branch of the “S” from “A'” to B (Fig. 0b) 
the importance of the orbits following the xl flow decreases. In 
Fig. [6b we can see the relative importance of the island around the 
x2-like stable orbit, centred at (x, x) = (3.56, -0.1823) with respect 
to the island around the xl-like, centred at (x, x) = (1.05, -8.2), at 
Ey — -32000. The local maximum of the histogram at about Ey 
« -32000 reflects the increasing importance of the stability island 
around the x2-like stable orbits as we approach this Ey value. 

This analysis shows that the observed ring structure in our re¬ 
sponse model is due to the x2-like orbits that are introduced in the 
system in the energy range where we have the “S” feature in the 
characteristic in the bistable region of the central family of periodic 
orbits. Stable periodic orbits with x2-like orientation exist on the 
xl' characteristic also beyond “B'” in Fig. 0b (up to Ey ^ -30000 
and provide to the system the backbone of the ring. The presence of 
the “S” feature in the characteristic is also associated with the end 
of the contribution of orbits trapped around stable xl-like orbits 
aligned with the y-axis to the response morphology of the model. 
As we have seen in Fig. [6b their contribution to the ring region 
is minor. Thus, the bar of the model ends practically at the ring. 
In general up to the energy of “A” and “A'” (Ey « -34500) our 
response model is populated by orbits associated with the lower 
branch of the xl' characteristic. Then we have a jump from “A'” to 
“A” and the backbone of our model are the stable orbits along the 
upper branch of the characteristic. 

3.2 The bar 

Eet us focus now in the response bar and its internal structure. If we 
confine ourselves to the study of the region of the response model 
(Fig. 0} with r < 2.5 kpc, we practically select the region of the 
bar. The majority of the particles with r < 2.5 kpc are located in 
the bar, since the regions to its sides, inside the ring, are regions 
almost depleted of particles. The energy distribution in this region 
is given in Fig. cit We observe that practically we have particles 
with Ey < -30000 and that the contribution of particles with Ey 
> -35000, i.e. to the right of A, A' in Fig. 0b, is small. The mode 
of the distribution is at Ey = -40750. 

The main feature of the response bar is an “X” feature, dis¬ 
cernible in Fig. \T\ We try to understand the mechanism leading to 
its formation first by investigating the possible contribution of par¬ 
ticles with the energy of the mode in Fig. la. The (x, x) surface of 
section at this energy is given in Fig. [7J). The stability island that 
dominates in the x > 0 region is the island of xl', which in this en- 
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Figure 8. The initial conditions of the integrated orbits on the (x, x) surface of section at Ey =-40750. Small empty bullets correspond to orbits supporting a 
narrow X feature matching the dimensions of the response bar. Small filled bullets correspond to orbits with an overall X morphology but not matching the 
dimensions of the bar. The large black bullets indicate the initial conditions of the orbits depicted in Fig. [9] 


ergy is xl-like. It is surrounded by a chaotic sea that extends also to 
x < 0. At the left side of this chaotic sea we observe two more sta¬ 
bility islands belonging to two families of periodic orbits that have 
a 3:1 character and their characteristic joins the branches of the xl r 
characteristic at energies beyond the “S” in Fig. |2] These orbits do 
not play any particular role in our study, so we will not deal further 
with them. 

In order to find out the orbits that support the “X” we con¬ 
sider a grid of initial conditions with a step 0.09 in the x- and 20 in 
the y- direction on the surface of section of Fig. We integrated 
the orbits with initial conditions on the nodes of the grid for 30 pat¬ 
tern periods. The initial conditions of the integrated orbits that have 
been found to give orbits with morphologies relevant to our study 
are given with filled and empty bullets in Fig. [8] They are encoun¬ 
tered for x > 0.5. In order to demonstrate their morphologies, we 
present below some characteristic orbital shapes. The large black 
bullets in Fig.[8]give the initial conditions of five typical orbits that 
are described in Fig. [9] All of them have xo=20 km s -1 , while their 
x 0 initial condition increases from Fig. [9^ to Fig.|9f. The large bul¬ 
let inside the innermost invariant curve in Fig.[8](xo=0.31) is given 
in Fig. [9^ and Fig.|9j). In Fig. [9^ it is drawn on top of the response 
model, so that we can see its extent relative to the bar and the ring 
structure. Its morphology can be seen in another scale in Fig. [9j> 
and is similar to the morphology of a xl-like periodic orbit. For 
x 0 =0.58, we are still in the area of the stability island (Fig. [9]e). 
The orbit has a boxy character, while we observe the appearance 
of wings of a X feature being formed towards the apocentra of the 
orbit. On the central stability island of Fig. [8] but close to its last 
KAM curve, the hole of the orbit around (0,0) almost vanishes and 
we have a fully developed X feature in the interior of the orbit. This 
can be seen in Fig. [9ji, with x 0 =0.67. Surrounding the stability is¬ 
land we find a plethora of chaotic orbits, that integrated at least for 
10 pattern periods reinforce the appearance of the X. In most cases 


the X feature was clearly discernible even after integrating the or¬ 
bit for 30 pattern periods. A typical case is the one given in Fig. 

For this chaotic orbit we have x 0 =1.03. The dimensions of this and 
all other similar chaotic orbits correspond to the dimensions of the 
response bar we see in Fig. Q] Finally, we find that starting inte¬ 
grating chaotic orbits at the outer border of the chaotic sea, close 
to the curve of zero velocity, we find again orbits with an X feature 
embedded in them. However, these orbits are rounder than all other 
orbits we discussed up to now. A characteristic example can be seen 
in Fig. [9f (x 0 =1.75). Due to their stubby morphology they do not 
contribute to the building of the structure of the response bar, since 
they extend further to its sides. 

Us ing an a lgorith m similar to that described in 
IChatzopoulos et alJ (1201 ll) . we applied simple criteria to character¬ 
ize the orbits according to their morphology and the degree they 
support the X-shaped bar of the response model in Fig.Q] In Fig. [8] 
small bullets indicate only the orbits we find supporting an X 
feature in the bar region. The empty small bullets indicate narrow 
orbits that remain confined in the region of the bar (|x| 1.5). Such 

orbits are like those depicted in Fig. E^,d and e. With filled small 
bullets we indicate orbits similar to the one in Fig. |9f. They have 
an interior X-shaped structure, though not as sharp as in Fig. [9^, 
but their dimensions do not match the dimensions of the bar. In 
conclusion the X-shaped bar is supported mainly by particles in 
orbits with initial conditions corresponding to the open bullets in 
Fig. [U 

By inspection of Figs. I7l8l andl9lwe first realize that the ini¬ 
tial conditions supporting the response bar morphology within a 
certain time are distributed on the surface of section independently 
of the locations of stability islands or chaotic seas. In general the 
stability region on Poincare cross sections that is excluded from 
participating in the building of a bar is found around the retro¬ 
grade family x4 (e.g. IContonoulos & Grosb0llll989!) . However, in 
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Figure 9. A series of non-periodic orbits with E j =-40750, integrated for 10 pattern periods, that demonstrate the orbital morphologies encountered at this 
energy. Their initial conditions are indicated with large bullets in Fig. [8] All of them have x'o=20 km s -1 , while from (a) to (f) their xo are 0.31, 0.31, 0.58, 
0.67, 1.03 and 1.75 respectively. The scale of the orbits is the same for panels (b) to (f). 


studies of barred-spira l potentials estimated from near-infrared ob¬ 
serva tions of galaxies (IChatzopoulos et al.l201 lklTsigaridi & Patsisl 
l2013h . it has been realized that there are initial conditions on xl 
stability islands, which do not support a particular morphologi¬ 
cal feature of a bar. On the other hand, we have found that there 
are initial condit ions in the chaotic seas that reinforce a particular 
bar morphology (IPatsis et al.l 19971l2010l : IChatzopoulos et al.l201ll : 

iTsigaridi & Pat sis 2013). In the present case quasi-periodic orbits 
that support part of the wings of an X feature are encountered at 
the periphery of the stability island of the xl-like periodic orbit. As 
noticed in IPatsis! (120051) the morphology of the periodic orbit at the 
center of an island may differ from the morphology of the quasi- 
periodic orbits of the outer invariant curves. For example, in cases 
of elliptical xl-like periodic orbits, the quasi-periodic orbits at the 
edge of the stability island may support an ansae-type morphology. 
As we see in Fig.[9]the morphology of the innermost quasi-periodic 
orbits, which is similar to that of the periodic orbit itself, does not 
resemble that of an X. On the other hand the integrated for 10-30 
pattern periods chaotic orbits with initial conditions in the chaotic 
sea surrounding the stability island have a morphology matching 
that of the response bar (cf. Fig. Q] and Fig. [9^). Since morpholo¬ 
gies similar to that of Fig. [9^ are encountered by integrating the 
initial conditions corresponding to the open bullets in Fig. [8] we 
can compare the area these open bullets occupy with the density of 
the consequents on the surface of section in Fig.^. It becomes ap¬ 
parent that the open bullets symbols are located in the sticky region, 
i.e. the region with larger density, surrounding the central stability 
island in Fig. Gk As regards the study of the orbital morphology it 
is worth to underline that moving from the center of the stability is- 



x (kpc) 


Figure 10. The periodic orbit of multiplicity 3 that has around it a sticky 
zone with chaotic orbits that have for several tens of pattern periods a mor¬ 
phology similar to the orbit in Fig.[9f. 

land outwards and then crossing the surrounding sticky region, we 
have a smooth morphological transition from narrow non-periodic 
orbits, along the x-direction, to broader non-periodic orbits, which 
harbour an X feature. The energy Ej =-40750 is the mode of the 
distribution in Fig. [7^. However, a similar analysis can be done for 



















Slowly rotating bars 9 



-10 -5 0 5 10 


x (kpc) 

Figure 11. Periodic orbits overplotted on the response model in the range 
-34000 <Ey < -25000. These orbits are the backbone for the support of 
the observed spiral pattern and the ring pattern. Energies from inside to out¬ 
side: -34000,-33500,-33000-32500-32000-31000-30000-29000, 
-28000, -27500, -27000, -26500,-26000,-25500 and -25000. 


every Ey roughly in the interval -45000 <Ey < -35000 (practi¬ 
cally for energies before the “S” feature). The non-periodic orbits 
supporting the X are mainly orbits in the sticky zones surrounding 
stability islands of xl-like periodic orbits. 

The zone with the small filled bullets in Fig.[8]is a sticky zone 
as well, this time around a periodic orbit of multiplicity 3. Its three 
tiny stability islands can be better seen in Fig.[7J) at (x, x) « (1.8,0), 
(-0.2, -130) and (-0.4,110). In Fig.[l0|we give this periodic orbit, 
so that its morphological relationship with the chaotic orbits in the 
sticky zone around it (Fig.(9f) becomes apparent. 


3.3 The spiral pattern of the model 

The response spiral is built by quasi-periodic orbits trapped around 
elliptical, precessing backwardly with respect to the direction of ro¬ 
tation as energy increases, xl' periodic orbits with Ey ^ -28000. 
We can see a set of these periodic orbits together with those that 
form the outer part of the ring in Fig. \Tl\ The mechanism that re¬ 
inforces the (trailin g) spiral arms is a typical “precessing ellipses 
flow” mechanism (IPatsisI 12009). In the energy range -28000 <Ey 
< -25000 the xl' stability island occupies practically all the (x, x) 
surface of section for v > 0. The orbits that support the spiral arms 
are in this case quasi-periodic orbits aroun d xl' . The orbital dy¬ 
namics are similar to those described bv IContonoulos & Grosb0ll 
CUE Il988h for open non-barred spirals. Corotation is placed at 
10 kpc, therefore the arms show clear bifurcating branches close 
the the inner 4:1 resonance. Close to r=8 kpc the xl' periodic or¬ 
bits have cusps and a clear rhomboidal morphology. The bifurcat¬ 
ing branches are formed by the congestion of these orbits with more 
circular orbits further away from the resonance. 


4 DISCUSSION 

We presented a dynamical mechanism that leads to the formation 
of a bar with a characteristic X embedded in it. In our model we 
have simultaneously also the formation of a ring tangent to the ends 
of the bar. We summarise this result in Fig. [12] We have used the 
model with Q p - 11.5 km s -1 kpc -1 for presenting this mechanism, 
because for this pattern speed the X feature extends all over the bar. 
This facilitates the description of the X morphology. However, it is 
a feature that we found also in the bar of other response models 
with Q p £ 15 km s -1 kpc -1 , without reaching the end of the bars. 
We give in Fig. El two snapshots from response models with Q p 
= 13.5 km s -1 kpc -1 (a) and 15 km s -1 kpc -1 (b) respectively. The 
maximum distances from the x-axis within which the X extends in 
these two models are indicated with dashed line segments parallel 
to the x-axis. The details of the X become less discernible as Q p 
increases. We observe that the X shrinks with increasing Q p from 
11.5 km s -1 kpc -1 (Fig. [TJ to 13.5 km s -1 kpc -1 (Fig. IT3ti) and 
then to 15 km s -1 kpc -1 (Fig. [T3b). This reflects the amount of 
sticky chaotic orbits that have been populated and participate in the 
building of the response bars in each case. In the models of Fig. [13] 
at distances beyond the end of the X feature, the bar has again a 
backbone of typical xl orbits until its end. 

The characteristic of the central family appears folded over 
a AE j « 2700 range, with both xl- and x2-like flows coexisting. 
The curve folds, but does not break. A similar behaviour has been 
encountered also in the slowly rotating “Model A2” of a 3D Fer¬ 
rers bar in lSkokos et al . ( 2002b ). The potential we study here is not 
a pure bar, but a barred-spiral one, allowing both bar and spirals 
appearing in the stellar response. Due to the low Q p , corotation is 
beyond the end of the spiral arms. 

Up to the E j of “A” and “A'” our response model is populated 
by orbits associated with xl-like periodic orbits (quasi-periodic or¬ 
bits from their stability islands and chaotic orbits from the sticky 
zones around these islands). For Ey larger than the Ey of “A” and 
“A'” our model is populated by orbits associated with the upper 
branch of the “S”. For the response model it is as “A” and “A'” co¬ 
incide on the characteristic. The presence of the bistable bifurcation 
creates an effective ILR, since it creates a x2-flow. The bistable bi¬ 
furcation determines the orbits that constitute the backbone of the 
response model. 

The ring still exists for Q p = 13.5 km s -1 kpc -1 (Fig. [I3k ) 
as a pseudo-ring, or double-ring structure. A rounder inner part 
seems to be detached from a more elongated x2-like outer part. 
In our main model (Fig.Q]) both parts appear joi nt forming a uniqu e 
ring. For Q p - 15 km s -1 kpc -1 (Fig. fl3b and iTsigaridi & Patsisl 
120131) attached to the ends of the bar we have only two density en¬ 
hancements of “smile” and “frown” morphology. However, a kind 
of loosely defined pseudo-ring structure appear even in models with 
Q p = 25 km s -1 kpc -1 . The ring in our main model is formed by or¬ 
bits having energies in the middle of the xl' characteristic in a (Ey 
,x) diagram (Fig. [2) and is associated with the folding of this curve. 
It is formed in the region where the central family of periodic orbits 
has a bistable character. The ring is formed due to the introduction 
in the system of a x2 flow, which in fact acts as if we had a usual in¬ 
ner Lindblad resonance region. By surrounding the bar however, it 
should be characterised as an “inner ring”, usual l y associated with 
inner 4:1 or 6:1 r esonances (see e.g. iButall 19861: iBvrd et al.lll994 
IPatsis et ai1l2003ll . Because of their association with inner n:l res- 
onances with n^ 4, inner rings are usually orientated along the bar 
t Butal 1986). In our case the “inner” ring can be characterised as 
almost circular or rather elongated towards the minor axis of the 
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Figure 12. (a) The central region of Fig. [I] which shows the inner part of the ring and the X feature in the bar of the model, (b) Two orbits from the sticky 
region around the xl' stability island in the surface of section for Ey = -39000 (x, x)=(0.25,-80) and (x, x)=(l,40) overplotted on the central region of the 
model. The corresponding xl' periodic orbit is xl-like, while the chaotic orbits in this sticky region give the orbital support to the X-shaped feature of the bar. 



Figure 13. Snapshots from the response models rotating with Q p =13.5 km s 1 kpc 1 (a) and 15 km s 1 kpc 1 (b). The dashed line segments parallel to the 
x-axis indicate the height within which extends the X in the two models. 


bar, depending on the energy level of the orbits we consider as the 
limit between the x2-flow orbits and the orbits that build the spiral. 

The X in the bar is mainly a result of populating sticky orbits 
around the xl' stability islands in the energy range for which we 
have a backbone of xl-like periodic orbits in our model. In princi¬ 
ple this can happen independently of the pattern speed of a model. 
However, in slowly rotating bars it appears pronounced, as it oc¬ 
cupies a larger part of them. In the case of the particular model in 
this paper the ring is tangent to an “X-shaped” bar. Our model de¬ 
scribes a configuration, where ordered and chaotic motion coexist 
and contribute to the formation of a unique structure. 

We note, that despite the 2D character of the present models, 
due to the morphology of the orbits that build the X, we can ob¬ 
serve a characteristic increase in the dispersion of velocities as we 


cross the ring inwards (Fig. d In Fig. d the arrows indicate a 
flow around the center of the model, while in Fig. fl4b the orienta¬ 
tion of the orbits looks random on the galactic plane. This is to be 
expected, since the motion of the particles following X-supporting 
orbits is complicated. 

iButa et all (120071 . pg.40). lEaurikainen et all J201 lL l2014h give 
several examples of X features in barred-spiral galaxies that are 
clearly far from edge-on (e.g. NGC 7020, NGC 1527, IC 5240, 
NGC 4429). The mechanism we propose here, i.e. the population 
of sticky chaotic orbits around the stability islands of the main fam¬ 
ily of periodic orbits, suggests an explanation of this morphology. 
We stress again that the X feature is pronounced in slowly rotating 
models, but slow rotation does not give the explanation of the fea¬ 
ture. The X appears as soon as the bar is built by sticky chaotic or- 
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Figure 14. (a) The velocity field of the model outside the ring area indicates a flow around the center of the system, (b) The velocity field of the model in the 
X-shaped bar region indicates an increased dispersion of velocities. 


bits in the chaotic zone surrounding the stability islands of xl-like 
periodic orbits, as well as by quasi-periodic orbits at the outer parts 
of these islands. The potential in the present paper originates in the 
estimation of the gravitational field of a late-type barred galaxy. 
Nevertheless, the X feature is expected to appear also in models 
for early type barred galaxies, provided that the bar is populated by 
the kind of sticky chaotic orbits we discuss. We note however, that 
there are cases of nearly face-on galaxies combining an X featur e 
with a ring (e.g. IC 5240, iButa et al1l2007l : lLaurikainen et alll20l4) . 
Such galaxies combine the morphological features of the response 
model of the present study. 

Recently Pats is & Katsanikasl J2014h investigated the orbital 
dynamics of 3D Ferrers bars that lead to the formation of boxy fea¬ 
tures inside the bars in their face-on projections. Despite the fact 
that the bars in the two cases are approached by different modelling 
techniques (2D vs. 3D dynamics), the two studies agree on the char¬ 
acter of the orbits that build the boxy features. In both cases we deal 
with sticky chaotic orbits. Similarly with what we find here, also in 
the 3D models the pool of orbits that are used to build the inner 
boxes are sticky chaotic orbits around the periodic orbits that con¬ 
stitute the backbones of the bars. However, the 3D modelling allows 
us to investigate the connection between the X features appearing 
in the face-on vi ews and the X featur es appearing in boxy/peanuts 
edge-on profiles dPatsis & Kats anikas 20 14|). 


5 CONCLUSIONS 

The main conclusions of our study refer to the dynamics of slowly 
rotating barred-spiral models and are enumerated below: 

(i) The characteristic of the main family in our slowly rotating 
barred-spiral potential includes periodic orbits that are successively 
(as E j increases) xl-like, x2-like, as well as elliptical orbits that 
precess if plotted at different Ej . It is characterised by a fold¬ 
ing, which we called the “S” (Fig. 0k). It is a case of two com¬ 
bined saddle-node bifurcations, usually called “bistability”. Along 
“S” the periodic orbits change their orientation. Slow rotation pro¬ 
nounces this folding. Our response model is populated by xl-like 
orbits up to the Ej of “A” and “A'” in Fig. UK and for larger Ej by 


orbits associated with the upper branch of the “S”. For the model it 
is as if we jump from “A” to “A'” in the characteristic. 

(ii) The abrupt change of the orientation of the response orbits 
leads to the formation of a ring in the middle of the characteristic, 
at the “S” region. The ring does not extend along the bar, as is the 
usual orientation of inner rings, but it is rather inclined towards the 
x-axis. We have effectively a x2 region beyond the end of the bar. 
The orientation and morphology of the ring resembles more that of 
nuclear rings. 

(iii) Bars built mainly by sticky orbits around the stability is¬ 
lands of the central family of periodic orbits have a boxy structure 
that harbours an X feature in it. In the class of barred-spiral mod¬ 
els we consider in the present study, the slower the pattern speed 
of the model, the greater is the importance of the X-feature for the 
overall morphology of the bar, since it occupies a larger fraction 
of it. However, the dynamical mechanism behind it is encountered 
in response models within a larger range of pattern speeds, since 
it depends on the amount of sticky chaotic orbits that populate the 
bar. 
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